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Abstract 

The set of 63 real generalized Pauli matrices of three-qubits can be factored into two subsets 
■ of 35 symmetric and 28 antisymmetric elements. This splitting is shown to be completely 

embodied in the properties of the Fano plane; the elements of the former set being in a 
bijective correspondence with the 7 points, 7 lines and 21 flags, whereas those of the latter 
set having their counterparts in 28 anti-flags of the plane. This representation naturally 
extends to the one in terms of the split Cayley hexagon of order two. 63 points of the 
hexagon split into 9 orbits of 7 points (operators) each under the action of an automorphism 
of order 7. 63 lines of the hexagon carry three points each and represent the triples of 
operators such that the product of any two gives, up to a sign, the third one. Since this 
hexagon admits a full embedding in a projective 5-space over GF(2), the 35 symmetric 
operators are also found to answer to the points of a Klein quadric in such space. The 28 
antisymmetric matrices can be associated with the 28 vertices of the Coxeter graph, one of 
two distinguished subgraphs of the hexagon. The PSL2(7) subgroup of the automorphism 
group of the hexagon is discussed in detail and the Coxeter sub-geometry is found to be 
intricately related to the ^-symmetric black-hole entropy formula in string theory. It is 
also conjectured that the full geometry/symmetry of the hexagon should manifest itself in 
the corresponding black-hole solutions. Finally, an intriguing analogy with the case of Hopf 



■ sphere fibrations and a link with coding theory are briefly mentioned. 
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1 Introduction 

Entanglement and geometry are key concepts in contemporary theoretical physics. The 
first can be regarded as the "characteristic trait of quantum theory" , the second as the 
basic unifying agent in explaining the fundamental forces. A vital combination of these 
concepts shows up in two seemingly unrelated fields, quantum information and quantum 
gravity. The basic goal in quantum information is understanding entanglement in geometric 
terms, and in quantum gravity the understanding of the structure of the Bekenstein-Hawking 
entropy of certain stringy black hole solutions. Recently, in a series of papers a striking 
correspondence has been established between this two apparently separate fields p]-[7]. For 
the intriguing mathematical coincidences underlying this correspondence M. J. Duff and S. 
Ferrara coined the term "Black Hole Analogy" . The basic correspondence of the analogy 
is the one between the black hole entropy formula of certain stringy black hole solutions 
on one hand and entanglement measures for qubit and qutrit systems on the other. The 
archetypical example PQ of such a correspondence is the one concerning the 8-charge STU 
black hole of N = 2, D = 4 supergravity [9]. The entropy formula in this case is of the 
form 

7T 

S=-^/tXbc, (1) 
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where tabc is the so called three-tangle [10], an entanglement measure characterizing en- 
tangled three-qubit systems A, B and C. Later works [21 HI H] have revealed that the most 
general type of stringy black hole solutions with 56 charges occurring within the context 
of N = 8, D = 4 supergravity also display this correspondence. Here the £7(7) symmetric 
black hole entropy formula is 

S = —^/tabcdefg, (2) 

where tabcdefg is an entanglement measure characterizing the tripartite entanglement of 
seven qubits. This measure is related to the unique quartic Cartan invariant for the fun- 
damental representation of the exceptional group £7(7), and to the discrete geometry of 
the Fano plane. Truncations of the discrete U -duality group £V(Z) have been connected 
to truncations of the Fano plane to its points, lines and quadrangles [5]. This result can 
be regarded as the first indication that discrete geometric ideas could play a key role in 
a quantum-entanglement-based understanding of duality symmetries. As a latest develop- 
ment, a physical basis for the Black Hole Analogy in the special case of the STU model has 
recently been revealed. It has been shown that the qubits occurring in the analogy can be 
realized as wrapped intersecting D3-branes [7]. 

It has become obvious from the very beginning that the quantum information theoretic 
basis for the black hole correspondence rests on the very special geometric properties of 
few qubit systems [3J. Indeed, in the quantum entanglement literature it has become com- 
mon wisdom that the special geometry of one, two and three-qubit systems can naturally 
be related to Hopf fibrations related to the complex numbers, quaternions and octonions, 
respectively [TI]-[T5J. These approaches to few qubit systems studied the geometry of en- 
tanglement via looking at the geometry of the corresponding multi-qubit Hilbert space. An 
alternative approach to understanding entanglement in geometric terms is provided by the 
one of studying the geometry of the algebra of observables with special operators (e.g., 
CNOT operations) creating and destroying entanglement by acting on it. In particular, it 
proved to be of fundamental importance to look at the structure of the so-called Pauli group 
on few qubits. This group is of utmost importance in the field of fault tolerant quantum 
computation and stabilizer codes. The algebra of the generalized Pauli group of iV-qubit 
systems has also remarkable discrete geometric properties, being completely embodied in 
the so-called symplectic polar spaces of order two and rank iV |14j-[17j. 

Although the appearance of few qubit (and qutrit) entangled systems related to discrete 
geometric structures provided an additional insight into the structure of stringy black hole 
solutions, some basic points still remained obscure. In particular, the entangled state charac- 
terizing the tripartite entanglement of seven qubits is not the one occurring in conventional 
quantum information theory. This state is not an element of any subspace of the Hilbert 
space on seven qubits. Indeed, it has been shown that it is an element of a 56-dimensional 
subspace of seven qutrits [4]. Moreover, a special role of three-qubit systems repeatedly 
showing up in these systems hints at an alternative structure lurking behind the scene. 

Motivated by the desire to eliminate such shortcomings, in this paper we initiate an 
observable-based understanding of the Black Hole Analogy. Our basic object of scrutiny 
is the discrete geometry of the Pauli group of three-qubits. Although the corresponding 
geometry of two-qubits, which is that of the generalized quadrangle of order two [151 116j. 
has already been worked out in detail, a similar comprehensive analysis for three-qubits is 
still missing. The main body of the paper is, therefore, devoted to a detailed study of the 
geometry of the real operators of this group. The principal result of this part is an explicit 
demonstration of the fact that the 63 real operators of the Pauli group can be mapped 
bijectively to the points of the so-called split Cayley hexagon of order two. The last part of 
the paper is devoted to establishing an explicit correspondence between a sub-geometry of 
this hexagon and the £7(7) symmetric black hole entropy formula. The symmetry of this 
(Coxeter-graph-based) geometry turns out to be related to a PSL2(7) subgroup of the full 
[/-duality group. In order to relate our observable based approach to the seven-qubit Hilbert 
space based one El [5], we provide the basic dictionary relating these two approaches. Finally, 
we also conjecture that the geometry and symmetry of the Hexagon in its full generality 
should manifest itself at the level of such black hole solutions. Our hope is that this new 
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approach will initiate further investigations to fill in the missing gaps, and to explore the 
possible physical consequences. 

The paper is organized as follows. First, we shall demonstrate that the algebra of our 63 
generalized Pauli matrices is fully encoded in the properties of the Fano plane and its dual 
(Sec. 2). Then, by employing its automorphism of order 7 (Sec. 3), we shall illustrate in many 
details how this encoding naturally and straightforwardly translates into the properties and 
characters of the points/lines of this split Cayley hexagon and one of its representation inside 
a projective 5-space (Sec. 4). Next, a quite in-depth examination of PSL 2 {7), a subgroup of 
the full automorphism group of the hexagon, will follow (Sec. 5). Finally, the Coxeter graph, 
one of two most important sub-configurations living within our hexagon (the other being the 
incidence graph of the Fano plane — the Heawood graph), will be given close scrutiny and 
shown to be intricately linked with the ^-symmetric black-hole entropy formula in string 
theory (Sec. 6). The paper will be finished by a brief discussion of an intriguing analogy 
with Hopf sphere fibrations and an outline of a worth-exploring link with coding theory. 



2 Real operators on three-qubits and the Fano plane 

Let us define the following set of 2 x 2 matrices 

Notice that except for Y these are the 2x2 identity and the usual Pauli spin matrices. 
They are symmetric except for Y which is antisymmetric. We note that the set M. = 
{±7, ±X, ±Y, ±Z} consists of the real operators of the Pauli group acting on a single qubit. 

We consider three-qubit systems with Hilbert space H = C 2 <g> C 2 ® C 2 on which the 
Kronecker (tensor) product of three 2x2 matrices A ® B <g> C is defined in the usual way. 
Here we refer to A, B and C as the matrices acting on the first, second and third qubit, 
respectively. Let us now assume that A,B,C € M. and introduce the shorthand notation 
ABC for their Kronecker product; hence, for example, we have 



ZYX =Z®Y®X 



Y®X 






-Y® X 



/ X \ 

-X 

-X 

\ OX / 



(4) 



Notice that the 8x8 matrices containing an odd number of Y"'s are antisymmetric, whilst 
those comprising an even number of Y's are symmetric; thus, for example, ZYX is an 
antisymmetric matrix, but IYY and XZX are symmetric ones. 

Up to a sign, we have altogether 4 3 = 64 distinct 8x8 matrices. After omitting the 
identity III = I ® J ® I we are left with 63 matrices of which 28 are antisymmetric and 
35 symmetric. In the space of symmetric 8x8 matrices we consider the following two 
seven-element sets 



C = {IIX, IX I, IX X, XII, XIX, XXI, XXX} (5) 

and 

V = {II Z, IZI, IZZ, ZII, ZIZ, ZZI, ZZZ} (6) 

to have a special status. We will see soon that the elements of these sets can be associated 
with the points (V) and lines (£) of the projective plane of order two, the Fano plane. Here 
we simply observe that the elements of these sets are pairwise commuting. 

Our next step is to create a multiplication table by taking the elements of the set C to 
label the rows and those of the set V to label the columns of a 7 x 7 array whose entries are, 
up to a sign, the (ordinary matrix) products of the corresponding elements. After carrying 
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Table 1: The matrix products between the operators from the two distinguished sets. 



Lines/Points 


ZZI 


ZII 


zzz 


IZI 


IZZ 


ZIZ 


IIZ 


IIX 


zzx 


ZIX 


ZZY 


IZX 


IZY 


ZIY 


IIY 


IXX 


ZYX 


ZXX 


ZYY 


IYX 


IYY 


ZXY 


IXY 


XIX 


YZX 


YIX 


YZY 


XZX 


XZY 


YIY 


XIY 


XII 


YZI 


YII 


YZZ 


XZI 


XZZ 


YIZ 


XIZ 


XXX 


YYX 


YXX 


YYY 


XYX 


XYY 


YXY 


XXY 


IXI 


ZYI 


ZXI 


ZYZ 


IYI 


IYZ 


zxz 


IXZ 


XXI 


YYI 


YXI 


YYZ 


XYI 


XYZ 


YXZ 


XXZ 



out all the calculations, we notice that the 21 symmetric and the 28 antisymmetric matrices 
from this 7x7 array reveal an incidence structure reminiscent of the one of the Fano plane. 
More precisely, we find three pairwise commuting symmetric and four pairwise commuting 
antisymmetric combinations occurring in both each row and each column. Moreover, any 
two elements such that one belongs to the symmetric and the other to the antisymmetric 
set (lying in a particular row or column) are anti-commuting (not commuting). Next, if we 
look at the pattern of the location of symmetric combinations pertaining to different rows 
and columns we obtain the same structure we get by looking at the incidence structure of 
points and lines of the Fano plane and its dual. 

Let us explicitly see how this works. After a particular ordering of the elements of the 
sets V and C we can group the remaining 49 matrices in the manner as depicted in Table 
1. In the first row (IIX) we have pairwise commuting symmetric combinations in columns 
1,2, and 4 (i.e., in the columns labelled by the triples ZZI, ZII, and IZI), and pairwise 
commuting antisymmetric ones in columns 3,5,6, and 7 (that is, in ZZZ, IZZ, ZIZ, and 
IIZ — the corresponding entries being underlined). Looking at the other rows we readily 
notice a cyclic shift of this pattern: (124) -> (235) -> (346) -> (457) -> (561) -> (672) -> 
(713). It represents no difficulty to check that this pattern is embodied in the structure of 
the Fano plane, as shown in Fig. 1, left; here, the grey numbers (points) respectively black 
ones (lines) from 1 to 7 label the columns respectively rows, together with the respective 
matrices of the two special sets, of Table 1. Dually, in the first column (ZZI) we have 
pairwise commuting symmetric combinations in rows 5,7, and 1 (XXX, XXI, IIX), and 
pairwise commuting antisymmetric ones in rows 2, 3, 4, and 6 (IXX, XIX, XII, and IXI). 
Now the cyclically shifted pattern looks as follows (134) -> (245) -> (356) -> (467) -> (571) 
— * (612) — > (723), being reproduced by the geometry of the dual Fano plane (Fig. 1, right). 




Figure 1: The Fano plane and its dual as the fundamental building blocks of the three-qubit 
Pauli group algebra. 
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IZI 



XII 



IXX (2 



ZZZ (3 



ZII 




5) XXX 



6)ZIZ 



IXI 



ZZI 



XXI 



Figure 2: The algebra of the operators from the two distinguished sets in terms of the 
incidence graph of the Fano plane. 

This dual view of the operators' algebra can be cast into a single compact form by 
employing the incidence graph of the Fano plane, also known as the Heawood graph (Fig. 2). 
This graph consists of 14 vertices corresponding to the seven points (grey circles) and seven 
lines (black circles) of the Fano plane, with two vertices (of different shading) being adjacent 
if they correspond to a point and a line such that the point is contained in the line. 

From the above-outlined picture we observe that the symmetric combinations correspond 
to the flags (a line taken together with a point lying on it) and the antisymmetric ones to 
the antiflags (a line and a point not lying on that line) of the Fano plane. Taken together 
with the sets V and £ whose elements label the points and the lines, we see that all 63 
nontrivial matrices forming a subset of the generalized three-qubit Pauli group can com- 
pletely be described in geometrical terms related to the structure of the Fano plane. This 
is a very important finding because there exists a remarkable unique finite geometry whose 
properties can be fully derived from the structure of the Fano plane alone, viz. a thick gen- 
eralized polygon called split Cayley hexagon [TB]-[52]. In order to show that this geometry 
indeed mimics all the essential features of the commutation algebra of our generalized Pauli 
matrices, one more piece of three-qubit mathematics must be given a careful inspection and 
properly understood — automorphisms of order 7. 

3 Automorphism of order seven and its orbits 

To this end in view, let us denote by a the following permutation: (1234567); this is obviously 
of order seven. We would like to find an 8 x 8 representation of a which permutes cyclically 
the elements of the ordered sets V and C we used to label the rows and columns of Table 1. 
We have to find, in particular, an 8 x 8 matrix T>(a) with the property (D(a)) 7 — III that 
acts on the elements of V and C by conjugation in the following way 

V(a) : IIX i > IXX ^ XIX ^ XII ^ XXX h-> IXI ^ XXI ^ IIX m- . . . , (7) 
V{a) : ZZI i > ZII ^ ZZZ ^ IZI ^ IZZ ^ ZIZ ^ II Z ^ ZZI (8) 



where 



V( a y 1 IIXV{a) = IXX, V{ a y 1 IXXV(a) = XIX, 



(9) 
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Such a matrix is of the form 



23(a) 



I P Q \ 

Q P 

QX PX 

\PX QX / 



P = 



1 




Q 




1 



(10) 



Now let us look at the action of 23(a) on 21 symmetric matrices labelling the flags and 28 
antisymmetric matrices labelling the anti-flags of the Fano plane. We find exactly seven 
orbits of seven elements each. The aggregate of 21 symmetric matrices splits into three 
orbits. One of them is represented by the elements located on the main diagonal of Table 1, 



23(a) : ZZX ^ ZXX ^ YZY ^ XZI ^ XYY ^ ZXZ ^ XX Z ^ ZZX i 



(11) 



the other two being represented by the matrices sitting at the positions shifted rightward 
by one column and/or three columns from the main diagonal in a cyclic manner, namely 



23(a) : ZIX ^ ZYY ^ XZX ^ XZZ ^ YXY ^ IX Z ^ YYI ^ ZIX h-s- . . . , (12) 



and 



23(a) : IZX ^ IYY ^ YIY ^ XI Z ^ YYX ^ ZXI i-> YYZ ^ IZX (13) 



respectively. The set of 28 antisymmetric matrices factors into four orbits; the elements of 
all of them are located above the main diagonal shifted from it by two, four, five and six 
columns. Their explicit forms read 



23(a) : ZZY ^ IYX ^ XZY ^ YIZ ^ XXY ^ ZYI ^ YXI ^ ZZY i 



(14) 



23(a) : IZY h-> ZXY ^ XIY ^ YZI ^ YXX ^ ZYZ ^ XYI ^ IZY i 



(15) 



23(a) : ZIY ^ IXY ^ YZX ^ YII i-> YYY h-> IYI h-> XYZ h-> ZIY h-> 



(16) 



and 



23(a) : IIY h-> ZYX ^ YIX ^ YZZ ^ XYX ^ IYZ ^ YXZ ^ IIY (17) 



respectively. 

It is important to realize here that the action of 23(a) on the thrce-qubit Hilbert space 
changes, in general, the entanglement type (see, e.g., the SLOCC classes a particular state 
belongs to [23 ). This can be demonstrated by expressing 23(a) in terms of three-qubit CNOT 
operations, which are well known to create or destroy entanglement 24 . One possible way 
to express 23(a) in terms of CNOT operations is 



23(a) = (Ci2C2l)(Ci2C3i)C , 23(Ci2C , 3i). 

Here, the explicit forms of the three-qubit CNOT operations are 



C 



12 



// 0\ 

2 

I 

\0 I 0/ 



C21 — 



(I 0\ 

I 

0/0 

\0 / 0/ 



(18) 



(19) 



c- 



23 



// 0\ 
0X00 
0/0 
X) 



C31 



(20) 



(P Q Q\ 
P Q 
Q P 
\0 Q P) 

with the first index relating to the location of the control bit and the second one to the 
target. 



G 



The fourth orbit of antisymmetric matrices, Eq. (|17[) . has several noteworthy proper- 
ties. The seven operators of this orbit are pairwise anti-commuting. Since the square of 
each of these operators is —III, replacing Y by 02 = —iY we get the matrices forming a 
representation of a seven dimensional Clifford algebra, 

{iTi, iT 2 , iT 3 , ir 4 , iT 5 , iT 6 , iT 7 } = {IIY, ZYX, YIX, YZZ, XYX, IYZ, YXZ} (21) 

a • I /. 1 , ^,,1. 1 = 111, j,k = 1,2,..., 7. (22) 

It is then straightforward to check that the remaining three orbits of antisymmetric matrices 
can be generated as the commutators of the form ^[Tj-,rfc]. This shows that the remaining 
21 matrices, Eqs. (JX4J) — (JT6j) , form the irreducible spinor representation of an so(7) algebra. 
Hence, the totality of 28 antisymmetric matrices Sok = fT/. and Sjk = j[Tj,Tk], associ- 
ated with the anti-flags of the Fano plane, form the generators of one of the irreducible 
spinor representations of the so(8) algebra with the generators Sjk giving rise to the so(7) 
subalgebra. 

Finally, we will notice that the elements of our specific fourth orbit can be used to la- 
bel the points of an oriented dual Fano plane. That is, instead of the ordered set C = 
{II X, I XX, XIX, XII, XXX, IX I, XXI} labelling the rows of Table 1 and the points of 
Fig. 1, right, we can use the ordered set £ = {IIY, ZYX, YIX, YZZ, XYX, IYZ, YXZ}. 
Since, for example, the result of the multiplication (XYX) (YZZ) (ZYX) (reading this 
from left to right) gives +IZZ, we can endow the line labelled by IZZ with an arrow 
pointing from left to right. Due to the above-described Clifford property, all even permu- 
tations of this triple product result in +IZZ, while odd ones yield —IZZ. Proceeding 
this way, every line of the dual Fano plane can be given a unique orientation — as il- 
lustrated in Fig. 3. Making one step further and employing the following correspondence 
between the unit octonions and the elements of our special set, {ei, e-2, €3, e±, e&, eg, ej} <-> 
{XYX,YZZ,YIX,ZYX,IIY,YXZ,IYZ}, Fig. 3 is found to be identical with a simple 
mnemonic for the products of the unit octonions [35]. Note, however, that for the octonions 
we have, for example, e\e 2 = e^, but in our case (XYX)(YZZ) ^ ZYX. So, the set C can 
only be used to mimic the signs occurring in the octonionic multiplication table. 




Figure 3: The oriented dual Fano plane. 

At this point we have at hand all the necessary technicalities to reveal the smallest 
split Cayley hexagon lurking behind the core algebraic properties of the 63 nontrivial real 
operators of the Pauli group on three qubits. 
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4 Core geometry: the smallest split Cayley hexagon 



A finite generalized n-gon Q of order (s, t) is a point-line incidence geometry satisfying the 
following three axioms [T8"]-|22j: 

(Al) Every line contains s + 1 points and every point is contained in t + 1 lines. 

(A2) Q does not contain any ordinary fc-gons for 2 < k < n. 

(A3) Given two points, two lines, or a point and a line, there is at least one ordinary 

n-gon in Q that contains both objects. 
From the definition it is obvious that the dual of a generalized n-gon is also a generalized n- 
gon. If s — t we say that the n-gon is of order s. An ordinary n-gon is the unique generalized 
n-gon of order 1. A generalized polygon is called thick if every point is contained in at least 
three lines and every line contains at least three points, i. e., if s, t > 2. The Fano plane is 
the unique generalized triangle of order two. There exists a unique (self-dual) generalized 
quadrangle of order two; this object that we already spoke about in the introduction plays 
a crucial role in the description of the generalized Pauli group of two qubits [HI [16] . It can 
be easily checked that the Heawood graph (the incidence graph of the Fano plane, Fig. 2) 
is the generalized hexagon of order (1, 2). There are also two generalized hexagons of order 
two, the so-called split Cayley hexagon and its dual [H]-[22], [26]. Our main concern here 
is the former one, also known as the G*2(2) hexagon because the Chevalley group G*2(2) is 
basically its automorphism group and in what follows simply referred to as the Hexagon. It 
has 63 points and 63 lines; every point is contained in three lines and every line contains 
three points. One of its stunning pictures created by the method of "finite pottery" [TBI [TU] 
is reproduced in Fig. 4. 

Since the Pauli group on three qubits contains 63 nontrivial real operators and the 
Hexagon is endowed with 63 points/lines as well, one is immediately tempted to conjecture 
that the two objects have something to do with each other. This conjecture is indeed fairly 
justified by the facts that the geometry behind the commutation algebra of the complex 
three-qubit Pauli group is that of the so-called symplectic polar space of order two and 
rank three (TU [17] and that our Hexagon lives fully embedded in this space [27] . Given 
the seven-fold symmetry of the Hexagon under which its 63 points are organized precisely 
into 9 orbits with 7 elements each (see Fig. 4) and the fact that we have an automorphism 
D(a) of order seven acting on 63 operators of the Pauli group which gives rise to the same 
orbit structure, all we have to do is to find the algebraic counterpart of the basic incidence 
relations of the Hexagon and verify the remaining structural relations. 

In order to carry out this program, we first notice that the Hexagon contains a copy of 
the Heawood graph (the sub-configuration in Fig. 4 consisting of the yellow and brown circles 
and the corresponding dark grey lines — compare with Fig. 2). According to our labelling 
scheme in terms of the ordered sets V and C, we notice that for the set {ZZI, IIX, ZII, 
IXX, ZZZ, XIX, IZI, XII, IZZ, XXX, ZIZ, IXI, IIZ, XXI} of consecutive points 
of this graph (Fig. 2) the operators corresponding to neighboring points commute. More- 
over, according to Table 1, after multiplying these commuting operators we get symmetric 
combinations (represented in Fig. 4 by the pink and orange circles) which again commute 
with their parents. In this way we obtain 28 points of symmetric combinations grouped into 
14 lines, each containing three pairwise commuting operators such that the product of any 
two of them yields (up to a sign) the third one. It is therefore reasonable to regard this 
property as the basic one of incidence governing the structure of our representation of the 
Hexagon. One further sees that the 14 lines we have just described are of the point /line/flag 
type (e.g., the triple {XXI,ZZI,YYI}). These lines are called Hl-lines [20], or lines of 
Coxeter type [27|. We have altogether 21 such lines (represented in Fig. 4 by dark grey 
arcs), with the remaining 7 of them connecting points related to commuting operators of 
point-line type to the remaining 7 symmetric operators (flags); these latter operators are 
represented by red circles in Fig. 4. Having employed the basic relations of incidence (aka 
up-to-a-sign multiplications producing triples of pairwise commuting operators), we have 
thus managed to establish the correspondence between 35 points of the Hexagon and 35 
symmetric combinations on the one hand, and between 21 Hl-lines and 21 triples of sym- 
metric combinations on the other hand. As a consistency check, one readily sees that the 
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XZY 




Figure 4: A diagrammatic illustration of the structure of the split Cayley hexagon of order 
two (after [18]-[2Qj) and how this geometry grasps the core features of the algebra of the 
real three-qubit Pauli matrices by bijectively associating the latter with the points of the 
hexagon. The points themselves are represented by circles whose interior is colored in nine 
different ways reflecting the nine orbits of an automorphism of order seven; the lines are 
drawn as black, light and dark grey curves/joints. The remaining symbols and notation are 
explained in the text. 
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order seven automorphism T>(a) rotates the operators corresponding to the yellow, brown, 
pink, orange and red points in the clockwise direction. 

We still have to take care of the remaining 28 points (anti-flags) of the Hexagon that 
answer to antisymmetric combinations of operators, as well as for the remaining 42 H2-lines 
[2D] , also known as lines of Heawood type (57 Let us deal first with the lines. These 
lines are of the nag/anti-flag/anti-flag type and their explicit construction is as follows. 
Take a (Fano) line L and a (Fano) point p lying on it. This gives rise to the flag {p, L}. 
Consider now the other two points q and r lying on the line L, and the other two lines 
M, N passing through the point p. Clearly, the two pairs {q, M}, {r, N} and {q, N}, {r, M} 
are two pairs of anti-flags. Hence, the flag {p, L} gives rise to two H2-lines of the form 
{{p, L}, {q, M}, {r, N}} and {{p, L}, {q, N}, {r, M}}. Since we have 21 flags the number of 
such lines is 42. In our representation in terms of operators of the Pauli group, these lines 
are represented by pairwise commuting triples containing two antisymmetric operators and 
a symmetric one with the property that any pair from the triple produces the remaining one 
under an up-to-a-sign multiplication. Given the representation of the Fano plane (Fig. 1, 
left, and Fig. 2) with its points labelled by the elements of the set V and its lines by those 
from C, it is a rather straightforward task to obtain all of these 42 triples, and to check our 
basic relations of incidence. As an illustrative example, one takes the point p represented 
by the operator ZII and the line L represented by the one IX X (see Fig. 1, left). From 
Table 1 one reads that they define the flag {p, L} represented by the operator ZXX. The 
remaining two points on L are q and r represented by IZZ and ZZZ, respectively. The lines 
M and N going through p correspond, respectively, to the operators II X and I XI. Hence, 
the operator representatives of the two pairs of anti-flags {q, M}, {r, N} and {q, N}, {r, M} 
are IZY, ZYZ and IYZ, ZZY, respectively. In this way we have obtained two H2-lincs 
represented by the triples {ZXX, IZY, ZYZ}, and {ZXX, IYZ, ZZY}. The elements of 
these triples are pairwise commuting and the product of any two from a given triple indeed 
yields, up to a sign, the remaining one. Finally, we observe that exactly one half of these 
lines (represented by black arcs/segments in Fig. 4) enjoy the property of connecting two 
distinct points/operators from the same orbit. 

When it comes to the remaining points of the Hexagon, all what we need is to make use 
of the special footing of the seven element set of antisymmetric matrices given by Eq. (fTT|) . 
For if we look at Fig. 4 carefully, the points (related to Fano anti-flags) represented by black 
bullets are also found to be special because the H2-lines passing through them (represented 
by light grey curves/arcs in Fig. 4) encompass not only all 21 points related to flags (red, 
pink and orange circles) but also all the remaining 21 points related to anti-flags (dark blue, 
light blue and green circles). None of the remaining points related to anti- flags exhibit this 
property. It is, therefore, natural to associate the points represented by black bullets with 
the matrices from C (Eq. (fT7]) ). Given this fact, and taking into account the structure of 
H2-lines, the seven-fold symmetry of the Hexagon and our up-to-a-sign multiplication recipe, 
after some experimenting we finally arrive at the full bijection between our real generalized 
three-qubit Pauli matrices and the points of the Hexagon as explicitly shown in Fig. 4. 

To gain a fuller appreciation of the above-described geometrical picture of our operator 
algebra, we have to mention a very important representation of the Hexagon as a sub- 
geometry of PG(5,2), the five-dimensional projective space over the Galois field of two 
elements [HI EH [29]. In this representation 63 points of the Hexagon coincide with 63 
points of PG(5, 2), 63 lines of the Hexagon form a special subset out of the total of 651 lines 
of PG(5, 2), and the 35 points of the Hexagon associated with the symmetric combinations 
of the operators are nothing but the points of a Klein quadric in PG(5, 2). And as there exits 
a bijection between this special set of points of PG(5, 2) and the 35 lines of PG(3, 2), which 
completely encodes this projective three-space (see, e.g., [30]), properties of our symmetric 
combinations are thus seen to be also embodied in the configuration of the lines of this 
smallest projective space. 

We shall conclude this section with the following important observation. Let us pick up 

lr The names of the two types of lines are dictated by the fact that after removing the points of flag type 
from the Hexagon we are left with two connected graphs — the Heawood graph and the Coxeter graph (see 
below); the edges of the former/latter corresponding to the lines of Heawood/Coxeter type. 
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one point represented by a black bullet. This point is associated with a unique set of seven 
points (including the point itself) which are located on the three lines passing through the 
point in question. In terms of our bijection, we thus have seven (one per each black bullet) 
distinguished sets of operators, viz. 



Bi 


= {YZZ, XX Z, ZYI, YXX, IYY, ZIY, XZX}, 


(23) 




= {YIX, ZXZ, XXY, YZI, IZX, XYZ, ZYY}, 


(24) 


03 


= {ZYXjXYYjYlZjXlYjYYZ^lYljZlX}, 


(25) 


Bi 


= {IIY, XZI, XZY, ZXY, ZXI, YYY, YYI], 


(26) 


B 5 


= {YXZ, YZY, IYX, IZY, YYX, YII, IX Z}, 


(27) 


B 6 


= {IYZ, ZXX, ZZY, XYI, XIZ, YZX, YXY}, 


(28) 


B 7 


= {XYX, ZZX, YXI, ZYZ, YIY, IXY, XZZ}. 


(29) 



Let us add to this list the two sets V and L, renamed as 

6 8 = V, B 9 = C. (30) 

After lengthy, yet straightforward, calculations one can verify that each of them consists 
of pairwise commuting operators and when taken together they give rise to a maximum 
(d +1 = 8 + 1 = 9) set of mutually unbiased bases (MUBs) in this particular three-qubit 
Hilbcrt space. 

5 The subgroup PSL2(7) of the full automorphism group 
of the Hexagon 

The Hexagon has a high degree of symmetry. By a symmetry we mean a bijection sending 
points to points such that this transformation induces a bijection between the line sets 
as well. We have already encountered the seven-fold symmetry of the Hexagon related to 
the permutation a and its 8 x 8 matrix representation T>(a) of Eq. (10). It is known that 
the full symmetry group of the Hexagon is of order 12096, being basically isomorphic to 
the Chevalley group G2(2) [19l[2T]. In this paper, however, we are merely interested in 
an important subgroup of the full automorphism group of the Hexagon realized on three 
qubits. This subgroup, which is directly related to the one of the Fano plane, is Klein's group 
PSL 2 (7) of order 168. In the following section, a 8 x 8 representation of this subgroup on 
three qubits will play an important role in illustrating how the geometry of the Hexagon 
manifests itself in the physics of certain stringy black hole solutions. 

The presentation of PSL2(7) we have found convenient [3T] is that based on three gen- 
erators a, (3 and 7 

PSL 2 (7) = {a, (3, 7 I a 7 = f3 3 = 7 2 = oT 2 ^- 1 = {iPf = (7«) 3 = !}■ (31) 

In the following two possible ways of realizing PSL 2 (7) in terms of permutations of the 
symbols 1, 2, . . . , 7 will be needed. They are 

a =(1234567), = (124) (365) (7), 7 = (12)(36)(4)(5)(7), (32) 

and 

a' = (1234567), = (1) (235) (476), 7' = (1)(3)(4)(25)(67). (33) 
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Our convention for multiplying permutations is from right to left (i.e., opposite to that 
adopted in Conway and Sloane [31 j ) . One can check that the generators a, (3, 7 and a' , (3' , 7' 
satisfy the defining relations of Eq. (|31|) . An important observation is that if we permute 
the points of the Fano plane by PSL2(7) transformations generated by the set a,(3,j then 
this implies a permutation on the lines generated by the elements a',/?', 7'. Hence, the set 
defined by Eq. ((32)) generates transformations for the points of the Fano plane and the set 
of Eq. (|33|) for the points of its dual (see Fig. 1). 

What we need is an 8 x 8 representation of PSL2(7) acting on three-qubits. It is clear 
that the 8x8 dimensional representation we are looking for is the one for which a is 
represented by 2?(a) of Eq. (10). For the remaining generators, we have found the following 
representatives 



T>Q3) = C 12 C 21 = 



(I 0\ 

/ 
0/00 

\0 / 0/ 



(34) 



D( 7 ) = C 2 i(/(8)/OZ) 



( z 








°\ 











z 








z 





V° 


z 





V 



(35) 



with the CNOT operations defined by Eq. (19). It is straightforward to check that T>(a), T>((3) 
and T>(-f) all satisfy the defining relations of Eq. (j3Tj) . It is important to realize that in this 
way we managed to realize Klein's group PSL 2 (7) on three qubits in terms of CNOT op- 
erations and a single "phase gate" Z. (The "phase" change in this case is just a sign 
change in the third qubit.) Hence, Klein's group operates quite naturally on the set 
M = {±/, ±X, ±Y. ±Z} of the real operators of the Pauli group. 

Since the order seven automorphism T>{a) is by its very construction a symmetry, in 
order to demonstrate the PSL 2 (7) symmetry of the Hexagon we have to check the action of 
the remaining generators T>{j3) and I? (7) on the three-qubit operators associated with the 
points of the Hexagon. To this end, it is convenient to replace our original labelling scheme 
(which was motivated by and based on the Polster and his co-workers' labelling of the points 
and lines of the Hexagon [18]-[20]) by the one shown in Table 2. So, for example, the set 

Table 2: The matrix products between the generalized Pauli operators from the two distin- 
guished sets in terms of the new labelling scheme (compare with Table 1). 



L/P 


hi 


h 2 


h 3 


hi 


h 5 


h 6 


h 7 


H 


ai 


bi 


Cl 


di 


ei 


fi 


9i 




92 


a-2 


b 2 


C2 


d 2 


e 2 


h 




h 


93 


a>3 


&3 


C3 


d 3 


e3 






h 


.94 


0,4 


iH 


C4 


di 








h 


95 


a 5 


&5 


C5 






d l:i 




h 


96 


Of; 


be 




h 


C7 


d 7 




h 


97 


a 7 



{gk I k — 1, 2, . . . , 7} corresponds to our set £' of operators of Eq. (17) playing a special 
role (black bullets in Fig. 4) and the lines of the Hexagon can be written in a nice compact 
notation as follows 



Ak — {fffc+ij c fc+3, a/c+4}, Bk — {fk-2, bk, 9k+i}, Dk = {dfe_i,pfe+i, efe+2}, 



(36) 



Cfe = {cfc_2, bk, Cfe+i}, Ek — {efc+i, efc+3, afc+4}, Fk = {dk-i, fk+x, /fc+2}, 



(37) 
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G k — {dk-i, hk+2, ik-i}, Hk — {ak+4,hk+4,ik+4},Ik — {bk,hk+i,ik}- (38) 

Just to familiarize ourselves with this new notation, we pick up one line of the Hexagon, 
say C\ = {cq, bx, C2}. According to Tables 1-2, this line is an H2-line of the form {ZYI, ZIX, 
IYX} depicted in Fig. 4 by the colour combination green-orange-green. In a similar way 
one finds that the line G 3 = {d,2,h^,i2} = {IYY,IZZ,IXX}, with the colour combination 
red-yellow-brown, is an Hl-linc. Obviously, the 42 lines of type A, B, C, D, E, F are H2-lines, 
and the 21 ones of G, H, I type are Hl-lines. The lines containing the special set C arc the 
ones of A, B, D; As already remarked in the previous section, (operators located on) these 
lines (represented in Fig. 4 by light grey acrs) give rise to a maximum set of MUBS. 

Let us demonstrate the PSL<2(7) symmetry of the Hexagon. Since we have already 
proved its seven-fold symmetry, we only have to check how the generators 2?(/3) and 2^(7) 
act on the points and lines of the Hexagon. Let us consider first the action of T>((3). Wc 
introduce the notation 

{rst) ijk = (r i s j t k )(sit J r k )(t t r j s k ), (rst)i = (nSiU), k = 1, 2, . . . , 7, (39) 

where (rst), as usual, refers to the cyclic permutation of the symbols rst. More generally, 
we can consider expressions like 

(rst)(uvw)( ijkKlmn)(p) = (rst) ljk (rst) lmn (rst) p (uvw) ijk (uvw) lmn (uvw) p (40) 

where (ijk)(lmn)(p) is a permutation of the ordered set {1, 2, . . . , 7} in a cycle notation. 
Terms like {r)ij k have no meaning hence they are not included. Thus, for example, (uvw)(r) ujiaq) 
stands for (uvw)ij k (uvw)iri. The action of T>(/3) on the points and lines of the Hexagon is 
thus given by 

(ad6)(ce/)(5)(»)(i)(253)(467)j (h) (142) (356) (7), (41) 

and 

(ADB) (CEF) (G/ff ) ( i42)(386)(7) , (42) 

respectively. Notice that this action is tied to the inverse of the basic permutations (3 and 
/?' of Eqs. (f3"2"|) (|3"3"|) related to the corresponding action on the sets V ({h k | k = 1, 2, . . . , 7}) 
and C ({i k | k — 1,2, ... , 7}) labelling the points and lines of the Fano plane. 

Let us now look at the action of T>{p/). A straightforward calculation shows that the 



points transform in the following way 

(ai&iXdxXci/iXeiJGfc), (a 3 d 3 )(b 3 )(c 3 )(e 3 )(f 3 g 3 ), (43) 

{a2d b ){d 2 a h ){b2b b ){c2gz){g2e- > ){e2fb){h c h), (44) 

(a 4 )(b4)(d 4 )(c 4 g4)(e 4 f4), (a 6 d 7 )(d 6 b 7 )(b e a 7 )(c e c 7 )(e e g 7 )(g e f 7 )(f 6 e 7 ), (45) 

{hih 2 ){h 3 h 6 )(h 4 )(h 5 )(h 7 ), (iiX^sX^X^X^?)- (46) 



It is important to stress here that the quantities in Eqs. (I43|) - (|44p . as well as i\, 12, is, and £5, 
transform with a sign change; thus, for example, (ai&i) and refer to the transformations 
a\ <-* —61 and i\ <-> However, since our labelling of the points of the Hexagon by the real 
Pauli operators is only up to a sign, this additional action of the group Z2 is not of current 
interest for us. Notice that the transformation rules given by Eqs. (|4"3"1) - (|4"5|) . as well as of 
those of the points {i k }, can be expressed by the permutations of the form (1)(3)(25)(4)(67), 
which is just 7' of Eq. (f33|) . Moreover, the transformation property (ab)(d)(cf)(e)(g) of 
Eq. (|4"3")l and the one of the points {h k } is associated with (12)(4)(36)(5)(7), i.e., with the 
permutation 7 of Eq. (f3"2"| . Of course, the transformations on the sets {i k } and {h k } exe 
again the ordinary ones of the points and lines of the Fano plane. 
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YYY 




Figure 5: The Coxeter graph, in a form showing its automorphism of order seven, as a 
subgraph/subgeometry of the Hexagon (see Fig. 4). 

Finally, we verify that the action of 2? (7) on the points also induces a bijection on the 



lines, which is of the form 

(A 1 D 3 )(A 2 F 1 )(A 3 C 6 )(A 4 C 1 )(A 5 F 6 )(A 6 D 4 )(A 7 ), (47) 

(B 1 S4)(S 2 B5)(B 3 C3)(S 4 C4)(S 6J B 3 )(B 7 C 7 )(C 2 C 5 )(C 7 F 7 ), (48) 

(D 1 E 2 )(D 2 F 2 )(D 5 F 5 )(D 6 E 5 )(E 1 F 3 )(E e F 4 )(E 7 ), (49) 

(G 1 H 2 )(G 2 )(G 3 i/ 1 )(G 4 i? 6 )(G 5 )(G 6 i/ 5 )(G 7 / 7 ), (50) 

(H 3 I 6 )(H 4 h)(H 7 )(I 2 I 5 )(I 3 )(h). (51) 



We see that under both T>((3) and 2^(7) the sets of HI- and H2-lines are left invariant. 
Hence, the PSL 2 (7) subgroup of the full automorphism group of the Hexagon preserves the 
splitting of its lines into these two subclasses. The same is also true for the two distinguished 
subsets of points. This is a very important finding which tells us that this subgroup also 
preserves the sub-configuration of the Hexagon formed by its 28 anti-flag points and 42 
H2-lines which is isomorphic — as depicted in Fig. 5 — to the famous Coxeter graph [32] , 
the second of the two remarkable graphs living in the Hexagon (see also [2D])- We shall 
make use of this result in the next section. For the sake of completeness, we here only add 
that the Hexagon contains 36 different copies of the Coxeter graph and the same number of 
copies of the Heawood graph as well. 

6 Coxeter sub-geometry and E^^-symmetric black hole 
solutions 

In this section we would like to show how the Coxeter sub-geometry of the Hexagon manifests 
itself in the geometry of a certain class of stringy black hole solutions. In string theory it 
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is well known that the most general class of black hole solutions in N = 8 supergravity/M- 
theory in four dimensions is defined by 56 charges (28 electric and 28 magnetic), and the 
entropy formula for such solutions is related to the square root of the quartic Cartan- 
Cremmer- Julia £7(7) invariant 1331 - 1361 



5 = 7TA 



Here, the Cremmer- Julia form of this invariant |36j 
8x8 central charge matrix Z, 

J 4 = Tr(ZZ) 2 - i(TrZZ) 2 + 4(PfZ + PfZ), 



(52) 

depends on the antisymmetric complex 

(53) 



where the overbars refer to complex conjugation. The definition of the Pfaffian is 
1 



PfZ 



ABCDEFGH 



ZabZcdZefZgh- 



2 4 ■ 4! 

An alternative form of this invariant is the one given by Cartan [34] 



■h = -Tr(a;y) 2 



4(Pfx + Pfy). 



(54) 



(55) 



Here, the 8x8 matrices x and y are antisymmetric ones containing 28 electric and 28 
magnetic charges which are integers due to quantization. These charges arc related to some 
numbers of membranes wrapping around the extra dimensions these objects live in 35j. The 
relation between the Cremmer- Julia and Cartan forms has been established [37j by using 
the relation 



1 



JJ 



4V2 



iyu)(r IJ ) AB , 



(56) 



where summation through the indices /, J is implied. Here (T ij )ab are the generators of 
the 5*0(8) algebra, where (I J) are the vector indices (I, J = 0, 1, . . . , 7) and (AB) are the 
spinor ones (^4, B = 0, 1, . . . , 7). Triality of 50(8) ensures that we can transform between 
its vector and spinor representations. A consequence of this is that we can also invert the 
relation of Eq. (|56|) and express x IJ + iyjj in terms of the matrix of central charge Zab ■ 

Since the matrices x IJ + iyu and Zab are both 8x8 antisymmetric ones with 28 inde- 
pendent components, it is natural to expect that the Coxeter sub-geometry of the Hexagon 
should provide additional insight into the structure of the £7(7) symmetric black hole en- 
tropy formula. In order to see that this is indeed the case, all we have to do is to recall from 
Sec. 3 that the elements of our distinguished set £ of Eq. (17), corresponding to the points 
of the Hexagon represented by black bullets (see Fig. 4), are related to a seven dimensional 
Clifford algebra. Moreover, commutators of the form [r^Tfe], j, k = 1, 2, . . . , 7, give rise to 
the 21 generators of an SO(7) subalgebra. Hence, the 28 antisymmetric matrices (T ij )ab 
occurring in Eq. (56) should be related to the 28 ones of the Coxeter set (see Eqs. (14)-(17) 
and Fig. 5). A natural choice for relating these two sets is 



{9k} 



(57) 



According to our labelling scheme of Table 2, we see that the 28 50(8) generators (T ij )ab 
are divided into the four sets {gk}, {bk}, {e/c} and {/&}, k = 1, 2, . . . , 7. An explicit relation 
between these quantities is 



( 


3i 


.92 


.93 


.91 


.95 


3o 


37 \ 


-.91 





e6 


c 4 


-h 


h 


-c 2 


-es 


-.92 







e7 


C5 


-h 


h 


-C3 


-33 


-c 4 







ei 


C6 


-h 


h 


-34 


fa 


-cs 


-ex 





e-i 


C7 


-h 


-35 


-h 


h 


-C6 


-e 2 





e3 


Cl 


-.96 




-h 


h 


-C7 


-ea 





e 4 


\-.97 




C3 


-h 


h 


-ci 


-e 4 


/ 



(58) 
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Here, the first row and column corresponding to the entries F ofe and r fc0 are related to the 
special set {gk}- The entries of the remaining 7x7 block are given by the formula h\j9j;9k]- 
Thus, for example, according to Tables 1 and 2 we have — T 25 = 5 [52 , 35] = \ [ZYX, XYX] = 
—YII = — / 4 . (All the entries of this matrix r /J are 8x8 matrices, hence the spinor indices 
AB are left implicit.) 

An important consequence of this explicit correspondence is that the 28 Gaussian integers 
(regarded as expansion coefficients for Z) comprising the antisymmetric charge matrix x IJ + 
iyu can also divided into seven sets. This division answers to the similar one of the "basis" 
vectors {gk, ek,Ck, fk}, k = 1, 2, . . . , 7, corresponding to T IJ . For historical reasons (to be 
explained below) we chose the dual labelling for the elements of the charge matrix x IJ + iyu 



J J 



yu 



Here, we used boldface letters for the coefficients ai 



/ 




-b 7 




-d 7 






&7 1 


a 7 





h 


d 4 


-c 2 


g2 


-b 4 


-ei 


b 7 


-fi 





gi 


e 4 


d 2 


a 2 


-c 4 


C7 


d 4 


-gi 





ai 


f 4 


-e 2 


b 2 


d 7 


c 2 


-e 4 


ai 





bi 


g4 


-fa 


e 7 


-g2 


d 2 


-h 


-bi 





Cl 


a 4 


f? 


b 4 


-a 2 


e-i 


-g4 


-ci 





di 




ei 


c 4 


b 2 


f 2 


a 4 


-dj 


/ 


fo 


-ao 


-b 


-c 


d 


-e 


-fo 


-go\ 


a 





f 6 


d 3 


-c 5 


g5 


-b 3 


-e 6 


bo 


-h 





g6 


e 3 


-d 5 


a 5 


-c 3 


c 


-d 3 


-g6 





a 6 


f 3 


-e 5 


b 5 


do 


C5 


-e 3 


-a 6 





b 6 


g3 


-f 5 


eo 


-g5 


d 5 


-f 3 


-b 6 





C<> 


a 3 


fo 


b 3 


a 5 


e 5 


-g3 


-c 6 





d 6 


\go 


e<> 


c 3 


-b s 


f 5 


-a 3 


-d 6 


/ 



(59) 



(60) 



, gj,J = 0,1,..., 7, not to be 



confused with the letters dj 
can be summarized as 



,gj, j = 1,2, ... ,7, denoting 8x8 matrices. This labelling 



{l,2,3,4,5,6,7}^{a,b,c,d,e,f,g} 



(61) 



{g,e,,c,f}~ {70,16,43,25}. (62) 

The symbols on the left-hand side of Eqs. (f6"Tj) - (|6"2")) label the basis vectors, those on the 
right hand side the expansion coefficients of the central charge matrix Z. In other words, 
we are converting numbers occurring in the sets {gk, ek,Ck, fk}, k — 1, . . . , 7, to boldface 
letters in a lexicographic order and the letters g,e,c, f to numbers according to the rule of 
Eq. (|62[) . This means that, for example, we use the notation for the expansion coefficient 
a; 36 + iy 3 6 attaching to the "basis vector" — T 36 ~ — f$ the symbol — e2 — ie$; that is, 5 was 
converted to e and the letter / was converted to the numbers 2 and 5 corresponding to the 
real and imaginary parts of the charge matrix. 

In this way we can associate the 28 points of the Coxeter sub-geometry of the Hexagon 
to 28 Gaussian integers with the real and imaginary parts corresponding to electric and 
magnetic charges characterizing the particular charge configuration of the black hole solu- 
tion. Explicitly, according to Tables 1 and 2, the sets of points represented by the (black, 
green, light blue, dark blue) circles in Fig. 5 can be associated with the "basis vectors" 
(<7fc, Cfc, efc, fk),k = 1,2,..., 7. Moreover, according to Eqs. (|5"8 j) - (|6"01) , with these "basis 
vectors" we can associate the corresponding "expansion coefficients," i.e., the charge com- 
binations. 

By virtue of the seven-fold symmetry of the Coxeter graph we also have managed to 
divide the 56 charges characterizing the particular black hole solution into 7 subsets with 
8 elements each. It can be shown [H [5] that each of these 7 sets can be regarded as 7 
three-qubit states with the 8 elements being their integer- valued amplitudes. The explicit 
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dictionary in the form of Eqs. ([59]) and (|60j) between the amplitudes of these three-qubit 
states and the elements of the charge matrices x IJ and yu were first obtained by Duff 
and his co-workers [38j . Here, we recovered this result by a straightforward way from the 
geometry of the Hexagon. It was also shown [U [5] that the transformation properties of the 
56 charges related to such three-qubit states under the fundamental representation of the 
exceptional group E 7 can be described by introducing an entangled state of seven qubits 
related to the geometry of the Fano plane. This seven-qubit system, however, contains 
merely tripartite entanglement in the form of our seven three-qubit systems we have just 
found. 

Notice that our present approach to this unusual type of entanglement is fundamentally 
different from those based on seven qubits. Here we managed to describe the charge con- 
figurations of black hole solutions by using the algebra of observables on i/iree-qubits. By 
comparison in these recent approaches such configurations have been described by special 
states on seuen-qubits forming a 56 dimensional subspace of the state space of seven qutrits 

®- 

In the previous section we have shown that the Coxeter sub-geometry has a PSL2(7) 
symmetry. Due to the relationship found between this geometry and the charge configu- 
rations we expect that this symmetry should also manifest itself in the black hole entropy 
formula. That this should be the case is already known from the mathematical literature 
|39j . However, now we can demonstrate this symmetry quite easily. Recall our 8x8 matrix 
representation of P5L 2 (7) generated by the matrices r D(a),'D(f3),'D('y). It is easy to check 
that these generators are elements of the group 5*0(8). Since under an element O € 5*0(8) 
the Pfaffian transforms as Pf (OZO*) = (DetO)Pf(-Z) = Pf(Z), from Eq. one can im- 
mediately check that J4, and so the black hole entropy formula, is invariant under P5L 2 (7). 
This should not come as a surprise since each of the terms of Eq. (|53|) has manifest SU (8) 
symmetry and all of our generators can be embedded into this subgroup of £7(7) ■ What 
is not trivial, however, is the explicit form of this P5L 2 (7) action on the charge matrix 
x IJ + iyu which is simply the one induced by a similar action of conjugation on the basis 
vectors T /J . This explicit action on the "basis vectors" {gt, Ck,ek, fk}, k = 1, . . . , 7, related 
to T /J via Eq. (58), is 

V{a) : (c)(e)(/)( 5 ) (1234567 ), (63) 



V{f3) : (ce/)(flO(i )( 253)(457), (64) 
T> (7) : (ci/i)(ei)(ffi), (c 2 35)(52e 5 )(e 2 /5)(/2C5), (c 3 )(e 3 )(/3.g3), (65) 

25 (7) : ( c 4.94)(e 4 /4), (c6C T )(e 6 g 7 )(g 6 f 7 )(f 6 e 7 ). (66) 

Here, we again note that the quantities of Eq. (|65ll ) transform with a sign change, i.e., 
( c i/i) means c\ «-> — /1 (see Eqs. (j4T |) . (j43^ -(f45" l) ). By virtue of the correspondence as given 
by Eqs. (|56p .([58 | - ((60|) . it is easy to obtain the corresponding explicit PSL^i^) action on the 
charge matrices x and y, or, alternatively, on the seven three-qubit systems built from seven 
qubits [U|5]. This explicit action of the P5L 2 (7) symmetry can be regarded as a subgroup 
of transformations of the full [/-duality group Ej(Z). Notice also that we have managed 
to represent this duality subgroup in terms of three-qubit CNOT operations creating and 
destroying entanglement and a trivial phase gate Z. 

As the last application of the Coxeter sub-geometry to stringy black hole solutions, let 
us consider the canonical form of the central charge matrix Z. As it is well known, by a 
transformation of the form Z 1— > U t ZU, where U € SU(8), Z can be brought into the form 



Z, 



canonical 



( Z1 








°\ 





Z\ 














Z3 













z&) 



Y, (67) 
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where all four Zj could be chosen to have the same phase, or three of the could be chosen 
to be real [35 . Observing that the antisymmetric combinations ci, ei, /i, gi occurring in the 
first row of Table 1 have the following explicit forms 



ci 



f 1 









-1 









-1 








V 



Y, ei 



(I 



Vo 





-1 









-v 



■ y, 



(68) 



/i = 



Z 1 







1 


\0 



\ 



-1 

-1/ 



Y, 



9i 



Y. 



(69) 



/l 0\ 
10 
10 
\0 1/ 

we see that obtaining the canonical form corresponds to the situation of attaching charges 
merely to the points of the Coxeter graph labelled by ZZY ', IZY , ZIY and IIY (see Fig. 5). 

Let us first assume that the numbers Zi are all real. Then, employing Eqs. (|56| . (|58p and 
US]), we get 



Vsz = x 01 gi + x 34 ei + x 57 Cl + x 26 h 
Hence, in this case we have 

- a 2 + ai + 84) 



-a 7 gi + aid + a 4 ci + a 2 /i. 



1 



Z1 = 7f (_a7 



z 2 = — =(-a 7 + a 2 - ai - a 4 ), 
V8 



1 



Z3 = —j=(-&7 - a 2 + ai - a 4 ), z 4 = — (— a 7 - a 2 - ai + a 4 ), 



V8 

with the entropy formula 



S = - 4aia 2 a 4 a 7 | 



(70) 



(71) 



(72) 



(73) 



being just the usual one obtained for four-charge extremal black holes [35]. (We remind that 
J 4 is positive for BPS and negative for non-BPS charge combinations.) 

In the usual "tripartite-entanglement-of-seven-qubits" interpretation [H [5J of the funda- 
mental representation of E7 we have seven three-qubit states with amplitudes a/, . . . ,gj, I = 
0, 1, . . . , 7. In this picture, obtaining the canonical form is simply restriction to just one of 
such three-qubit states (in this case to the one with amplitudes a/). This process is equiv- 
alent to the one with — a 7 — zao attached to the black circle labelled by IIY, and a 4 + 
&2 +i&5 and ai to the circles labelled by ZZY, ZIY and IZY, i. e., to the correspond- 
ing green, dark blue and light blue elements of Fig. 5. As explained elsewhere [H IH [5J, the 
black hole entropy formula in this case is related to Cayley's hyperdeterminant [I], i.e., to 
the unique triality and SX(2)® 3 invariant entanglement measure for three qubits. Indeed, a 
straightforward calculation in this case shows that J 4 = — D(a) where the explicit expression 
for Cayley's hyperdeterminant D(a) is 



D(a) = (a a 7 ) 2 + (aia 6 ) 2 + (a 2 a 5 ) 2 + (a 3 a 4 ) 2 -2(a a 7 )[(aia 6 ) 
- 2[(aia 6 )(a 2 a 5 ) + (a 2 a 5 )(a 3 a 4 ) + (a 3 a 4 )(aia 6 )] 
+ 4a a 3 a 5 a 6 + 4aia 2 a 4 a 7 



(a 2 a 5 ) 



(a 3 a4)] 
(74) 



giving a generalization of Eq. (J73J as S = — D(a)\. 

Obviously the action of 2?(q;) on the canonical form of Z by conjugation gives rise to 
the permutation of the charges in the form (abcdefg). Combining this with the P5L 2 (7) 
invariance of J 4 shows that in agreement with the well-known result [H [5J [3j|] the expression 
for J 4 in terms of the 56 charges contains (among others) the sum of seven copies of Cayley's 
hyperdeterminant. These terms describe the seven different possible N — 2, D — 4 STU 
truncations of the N — 8 solutions. Using the explicit dictionary of Eqs. (f58]i — (f60 |) we can 
alternatively label the points of Fig. 5 by the 28 Gaussian integers. It is easy to check that 
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the charge configurations belonging to a particular truncation of this type correspond to 
seven sets of four points with mutual distance four. 

After we have clarified the geometrical meaning of the canonical form for 2, let us try 
to understand its symmetry properties. As already explained, we have only eight charges 
corresponding to the eight integers aj, I = 0, 1, . . . , 7. The action of the PSL 2 (7) subgroup 
of the duality group on these charges can be understood from the corresponding action on 
the operators ci,ei,/i,5i (see Eqs. ([6"T j) - ([62|) ). From Eqs. (|64|) and ([65]) we find that V(j3) 
and I?(7), acting as (ciei/i)(gi) and (ci/i)(ei)(<7i), generate the dihedral group D3, which 
is isomorphic to the permutation group S3. The action of S3 on the eight charges in the 
permutation notation reads 

£>(/?) : (a )(a7)(a3a 6 a5)(aia 2 a4), (75) 

and 



X>( 7 ) : (a )(a 7 )(a 1 )(a 6 )(a 2 a4)(a 3 a5). (76) 

At this spoint, we get in touch with the usual interpretation based on seven three-qubit 
systems |H |SJ- Let us associate with the 8 charges aj, I = 0, 1 . . . , 7, the three-qubit state 

\i/>) =aooo|000) +800x1001) + ..- + a U i|lll), (77) 

i.e., we reinterpreted the charges a/ (decimal labeling) as integer-valued amplitudes a. a bc, 
a,b,c G {0, 1} (binary labeling) of a three-qubit state. Then the transformation rules of 
Eqs. (|75|) and (f76|) for the amplitudes correspond to the permutations (321) and (12) of 
subsystems (we are labeling qubits from left to right). Since these permutations generate 
the full permutation group on three letters, the corresponding truncation of J4 (i. e., Cayley's 
hyperdeterminant) should be a permutation invariant, and this is indeed the case. Hence, as 
another important consistency check, we reproduced the permutation invariance of Cayley's 
hyperdeterminant from the PSL 2 (7) symmetry of Cartan's quartic invariant. 



7 Conclusion 

We have performed a comprehensive analysis of the properties of the real generalized Pauli 
operators of three-qubits to the extent comparable with the treatment of the (complex) 
two-qubit case [El [16] . The algebra of the operators is found to be completely describable 
in terms of the structure of the Fano plane and uniquely extendible to the picture of the 
split Cayley hexagon of order two. The 63 Pauli operators are bijectively identified with 63 
points of the Hexagon whose 63 lines carry each triples of operators; these triples are such 
that the product of any two of them gives, up to a sign, the third one. A deeper insight 
into the structure of the Hexagon is acquired in terms of one of its automorphism of order 
7 that features 9 orbits of 7 points/operators each. The factorization of the set of operators 
into 35 symmetric and 28 antisymmetric members finds its natural explanation in terms 
of two distinct kinds of points of the Hexagon: 7 Fano points, 7 Fano lines and 21 Fano 
flags on the one hand versus 28 Fano anti-flags on the other hand, respectively. Triples of 
collinear operators are of two distinct kinds as well, according as they lie on the lines of 
point/line/flag (Heawood) type or on those of flag/anti-flag/anti-flag (Coxeter) type. The 
geometry behind 28 antisymmetric matrices is governed by the properties of the Coxeter 
graph and its symmetries. The geometry behind 21 out of 35 symmetric guys rests on the 
Heawood graph, while that of the full set is mimicked by a Klein quadric in a projective 
5-space over GF(2). 

Employing this novel finite geometrical language, an alternative quantum informational 
view of the structure of the £V-symmetric black hole entropy formula and a subgroup of its 
duality symmetry has been arrived at. This view is based on the algebra of observables on 
three-qubits. This is to be contrasted with the other ones [11 [5] based on an entangled state of 
seven qubits. In both pictures, however, three-qubits play a distinguished role. We managed 
to establish a dictionary between these complementary approaches. However, there are still 
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many questions left to be answered. The most exciting one is the possibility of relating the 
black hole entropy formula to the full geometry of the Hexagon. Indeed, we only managed to 
relate black hole solutions to the Coxeter sub-geometry with an explicit PSL2(7) symmetry. 
We conjecture that the full C?2(2) symmetry of the Hexagon should play a role in this black 
hole context. One solid piece of evidence in favour of this conjecture is the fact that since E? 
is of rank seven we have 133 — 7 = 2 x 63 = 126 step operators. According to an important 
result of M. Koca and his co-workers [40], the adjoint Chevalley group G2(2) acts naturally 
on the (octonionic) root system of E-j, Using this result and our explicit dictionary, it would 
be really spectacular to establish a full correspondence between the geometry of the Hexagon 
and the duality symmetries of stringy black hole solutions. 

As already remarked, there exist two distinct generalized hexagons of order two: our 
Hexagon and its dual. One is naturally tempted to ask why three-qubits favour the former, 
not the latter. Since the algebra of the corresponding operators is encoded in the structure of 
the Fano plane, the answer seems to be rather easy: our Hexagon, as we have seen, contains 
a sub-hexagon of order (1, 2) which is isomorphic to the incidence graph of the Fano plane, 
but its twin does not. Yet, there exists another remarkable structure within our Hexagon 
which is absent in its counterpart, a distance- 2- ovoid |18j- 22 . A distance- 2-ovoid is a set 
of non-collinear points such that every line is incident with exactly one point of that set 
and for a generalized hexagon of order (s, t) its cardinality amounts to s 2 t 2 + st + 1. In our 
Hexagon, this is the unique configuration of 21 points of flag type (orange, red and pink 
circles in Fig. 4) . Why we should be bothered with this set? Because generalized hexagons 
and their sub-configurations have recently been found to be an important source of a variety 
of codes [41j . with those arising from distance-2-ovoids being of specific two- weight type. 
This link with coding theory is something which certainly deserves a serious perusal of its 
own. 

As a final note, we would like to draw the reader's attention to the following puzzling 
analogy. We have already mentioned in the introduction that multi-qubit states of smallest 
orders can geometrically be modelled by Hopf fibrations. Namely, the fc-qubit states, k = 
1, 2, 3, are intimately connected with the Hopf sphere fibrations of the 

o2(* ! + 1 '-i s 2 "- 1 S 2 k 

type, respectively |11[ 113] . A very similar situation occurs in the setting provided by thick 
generalized n-gons. For if we associate the geometry of three non-trivial matrices given by 
Eq. (3) with that of the generalized digon (n = 2) of order two, and take into account that 
the two-qubit case is underlaid by the geometry of the generalized quadrangle (n = 4) of 
order two p~5l[16], we also arrive at such a three-fold correspondence: n = 2k, k — 1,2,3. 
Moreover, in both the cases, remarkably, the sequence ends at k = 3; in the former case 
due to the fact that there are no natural Hopf fibrations of higher order, whilst in the 
latter one because, by the celebrated result of Feit and Higman [42], of non-existence of 
thick generalized n-gons with n > 6 and meeting the constraint s — t. We thus observe a 
fascinating phenomenon when two substantially different pieces of mathematics convey the 
same message: a special footing played by single-qubit, two-qubit and three-qubit systems 
within an infinite family of multiple-qubit states. Is this analogy a mere coincidence, or is 
there a deep physical reason behind? 
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